Azimuthal dependence of pion source radii 
in Pb+Au collisions at 158 A GeV/c 
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We present results of a two-pion correlation analysis performed with the Au+Pb collision data 
collected by the upgraded CERES experiment in the fall of 2000. The analysis was done in bins 
of the reaction centrality and the pion azimuthal emission angle with respect to the reaction plane. 
The pion source, deduced from the data, is slightly elongated in the direction perpendicular to the 
reaction plane, similarly as was observed at the AGS and at RHIC. 

PACS numbers: 25.75.Gz, 25.75.Ld 



I. INTRODUCTION 

Two-particle correlations provide unique access to the 
spatial extension of the source of particles emitted in the 
course of heavy ion collisions at relativistic energy (for 
a recent review see The relation between the ex- 

perimental correlation function, defined as the relative 
momentum distribution of pairs, normalized to the anal- 
ogous distribution obtained via event mixing, and the 
size of the fireball is especially simple in case of identical 
pions where the main source of correlations is the Bose- 
Einstein statistics. In fact, in the pure boson case the 
correlation function is 



C(q,P) = l 



1/ d*x S{x, P) exp(iq ' 
\JdixS{x,P)\^ 



(1) 



with the source function S{x,P) describing the sin- 
gle particle density in 8-dimensional position-momentum 
space at freeze-out. The correlation function C depends 
on the momentum difference q = P2 - Pi and the mean 
pion momentum k = (pi+p2)/2. The width of the peak 
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at q = is inversely proportional to the source radius. 
A particularly exciting prospect is to look for a source 
asymmetry possibly refiecting the initial asymmetry of 
the fireball created in collisions with finite impact pa- 
rameter. Indeed, significant dependence of the source 
radii on the azimuthal emission angle with respect to the 
reaction plane, defined by the beam axis and the impact 
parameter vector, was observed in Au-|-Au collisions at 
2-6 GeV and at yG=200 GeV 0. In this paper we 
present results of the first analysis of the azimuthal de- 
pendence of pion source radii at SPS energies. 



II. EXPERIMENT 

The CERES/NA45 experiment at the CERN SPS is 
described in detail in [j]. The upgrade by a radial 
Time Projection Chamber (TPC) in 1997-1998, in addi- 
tion to improving the dilepton mass resolution, enhanced 
the experiment's capability of studying hadronic observ- 
ables. The cylindrical symmetry of the experiment was 
preserved during the upgrade, making the setup ideally 
suited to address azimuthal anisotropics. About 30 mil- 
lion Pb-}-Au collision events at 158 A GeV/c were col- 
lected in the fall of 2000, most of them with centrality 
within the top 7% of the geometrical cross section ctg- 
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Small samples of (t/(tg=20% and minimum bias colli- 
sions, as well as a short run at 80 A GeV/c, were recorded 
in addition. 



III. DATA ANALYSIS 

The results presented here are based on a correlation 
analysis of the high-statistics Pb-|-Au data set from the 
year 2000 The results of an azimuthal-angle aver- 
aged analysis in are consistent with the previously 
published CERES data and with a recent analysis 
by the NA49 Collaboration 0. The main steps of the 
azimuthal-angle dependent analysis are described below. 



calculated (modulo tt) from the calibrated Q2 via 

1 

$RP = ^arctan^. (4) 

The resolution of the so determined reaction plane angle, 
estimated via the subevent method [1, was 30-34°. 

The event mixing, needed to obtain the denominator 
of the correlation functions, was performed in bins of cen- 
trality (2% of ctq), event plane (15°), and vertex (same 
target disk). This ensures that the shape of the back- 
ground is identical to that of the signal in all respects 
except for the femtoscopic correlations. 



A. Event selection 

The collision ccntrality was determined via the charged 
particle multiplicity around midrapidity ?/boam/2=2.91. 
Two variables, the amplitude of the Multiplicity Counter 
(MC) (single scintillator covering 2.3< 77 <3.5) and the 
track multiplicity in the TPC (2.1< 77 <2.8), were al- 
ternatively used as the centrality measure. Knowing the 
data acquisition dead time factor and the target thick- 
ness, and assuming that all beam particles were hitting 
the target, the event counts were translated to the cross 
section for collisions with a given multiplicity. The cen- 
trality variable used in this paper is defined as the in- 
tegrated cross section divided by the geometrical cross 
section (Tg=6.94 barn. 

The fireball created in a collision with finite impact 
parameter is elongated in the direction perpendicular to 
the reaction plane. In the course of expansion, with the 
pressure gradient larger in-plane than out-of-plane, the 
initial asymmetry should get reduced or even reversed. 
Experimentally, the source eccentricity at the time of de- 
coupling can be determined from an analysis of two-pion 
correlations as function of $* = <&pair — ^rPi the pair 
emission angle with respect to the reaction plane. 

The azimuthal angle of the reaction plane <i>Rp was 
estimated by the preferred direction of the particle emis- 
sion (elliptic flow). For this, in each event the flow vec- 
tor Q2 was constructed out of the measured particles, 
characterized by transverse momenta pt and azimuthal 
emission angles (/) [1, Q : 

1 ^ 

= ]^IIpm cos(20O (2) 
i=i 

1 ^ 

= N^^'-' ^i^(20.)- (3) 

1=1 

The Q2 and Q\ components were calibrated (shifted and 
scaled) such that the peak in the {Q2 tQ^. ) distribution 
was centered at (0,0), and its widths in the directions 
of X and Y were equal. The reaction plane angle was 



B. Track selection 

Only TPC tracks with at least 12 (out of maximally 
20) hits and a reasonably good x^j falling into the fiducial 
acceptance 0.125 < 9 < 0.240, were used in the analysis. 
A momentum dependent dE/dx cut was applied to re- 
duce the contamination of the pion sample. Pions from 
and A° decays were suppressed by a 2.5 a matching 
cut between the silicon vertex detectors and the TPC. 



C. Pair selection 

The two-track resolution cuts applied to the true pairs 
and to the pairs from event mixing were different for 
the two possible track pair topologies: in the case of the 
magnetic field bringing the tracks apart from each other 
(sailor) the required two-track separation was A0 >38- 
45 mrad, depending on the transverse momentum; for 
the opposite (cowboy) case, Acj) >90-140 mrad was used. 
The polar angle cut was the same for both topologies 
(A6'>8-9 mrad). It should be noted that the required 
two-track cuts depended somewhat on the quality cuts 
applied to single tracks: the higher the number of hits 
required for single tracks, the more pairs were lost be- 
cause of the finite two-track resolution. 



D. Correlation functions 

The two-pion analysis was performed in the longitudi- 
nally co-moving frame (LCMS) defined by the vanishing 
z-component of the pair momentum. The momentum 
difference in this frame, q = P2-P1, was decomposed into 
the "out" , "side" , and "long" components following the 
Bertscli-Pratt convention, with gout pointing along the 
pair transverse momentum and giong along the beam [ll[ . 
The particles were numbered such that gside was always 
positive. 

The 7r~7r~ and tt+tt^ correlation functions were fitted 
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by 

C2(q) = TV • {(1 - A) + A • ^^4q) (1 + G(q))} , 
G(q) = exp(-^i?,2,.q,g,) (5) 

with the indices i,j = {out, side, long}. The normahza- 
tion factor was needed because the number of pairs 
from event mixing was four times higher than of signal 
pairs. The correlation strength A < 1 reflects the con- 
tribution of pions from long-lived resonances, the con- 
tamination of the pion sample by other particle species, 
and the finite q-resolution. The resulting source radii 

^jRij describe the size of the source emitting pions of 
a given momentum The Fc{qinv) factor, gi„„ = 

•\/—(j)2 ~ Pi accounts for the mutual Coulomb inter- 
action between the pions and was calculated by averaging 
the nonrelativistic Coulomb wave function squared over 
a realistic source size. The Coulomb factor was attenu- 
ated by A similarly as the rest of the correlation function 
peak; the importance of this approach was demonstrated 
in . The fits were performed by the minimum negative 
log-likelihood method, assuming that the number of true 
pairs in a given bin is distributed around the expected 
mean value (equal to the number of mixed pairs times 
the fit function) according to the Poisson statistics. The 
source radii obtained from the fit were corrected for the 
finite momentum resolution 



Ap 

P 



= 2%el%-p/(GeV/c). 



(6) 



The correction was determined by Monte Carlo and is 
rather insignificant for i?side and i?iong; for i?out it gets 
as large as w 4-20% for the highest pion momenta. 



E. Azimuthal dependence of pion source radii 

Correlation functions were generated separately for 
pion pairs with different azimuthal angles with respect 
to the reaction plane $* = $pair ^ ^rp- For this, the 
pions were sorted into eight bins covering (— 7r/2, tt/2). 
The eight correlation functions were fitted [1^| and the 
resulting squared source radii were parametrized by 



-"-out 
n2 

-'''side 
d2 

-"-long 

n2 
-'''OS 

■Rol 



^out.0 + 2i?out,2COs(2$*) 
i?2dc.0 + 2i?sido,2COs(2$*) 
^lo„g,o + 2Aong.2COs(2$*) 

2i?2 sin(cf>*). 



(7) 



While the Rl^i Q, Rfi^^ q, and Rf^ng o obtained represent 
the phi-averaged squared pion HBT radii, the second 
Fourier coefficients /Jfj describe the eccentricity of the 



observed pion source. Since the reaction plane is known 
modulo TT odd Fourier components should vanish. If the 
pion source were to reflect the initial collision geometry 
(almond shape out-of-plane) a positive R^^^g 2 -^os 2 
and a negative i?out should be expected. For symmetry 
reasons all anisotropics should disappear in the limit of 
central collisions. 

The second Fourier coefficients R^2 have been cor- 
rected for the reaction plane resolution by dividing them 
by the mean cosine of twice the difference between 
the reconstructed and the true reaction plane angles 
(cos[2(^'g'p — 5'rp)]), estimated via the measured mean 
cosine of the difference between two subevents: 



Rf-- 



v/2(cos[2(*f, ^ 



(8) 



similarly as it is done for flow measurements [T3|- The 
correction factor was between 4.8 and 2.4 (centralities 
of 0-2.5% and 25-70%, respectively). The appropriate- 
ness of the flow correction for two-pion correlation radii 
has been verifled with a Monte Carlo simulation of pion 
pairs emitted from an elliptical source; the observed re- 
duction of the second Fourier coefficients of squared radii 
was, within the statistical errors of the simulation (about 
10%), in agreement with Eq. ([5]). 



F. Systematic errors 

The fact that separate analyses of positive and nega- 
tive pions give consistent results indicates that the lim- 
ited particle identification is not causing any bias. The 
systematic error related to the reaction plane resolution 
correction, based on the comparison between Eq. ([H]) and 
a numerical simulation, is estimated to be not larger than 
10%. Autocorrelations, i.e. using the same particles to 
determine the event plane and the HBT radii, might be 
another source of systematic errors; unfortunately, an at- 
tempt to exclude the two particles from the event plane 
construction led to a strong distortion of the correlation 
functions as it was not clear how to do the exclusion con- 
sistently for the true and mixed pairs. An independent 
way to estimate the overall systematic errors is by in- 
specting (see Fig. [T] in the next section) the coefficients 
i?oi^i and i?si^i which should be zero because the reaction 
plane is known only modulo tt. Averaged over central- 
ities, they are 0.40 (16) and 0.22 (6) fm^, respectively. 
Attributing part (one sigma) of the deviation to statis- 
tical fluctuation we are left with a discrepancy of about 
0.3 fm^ which we thus assume to be the systematic un- 
certainty of the results presented. 



IV. RESULTS AND DISCUSSION 

The extracted Fourier coefficients are shown in Ta- 
ble |T] and Fig. [1] [1^. The anisotropics in the out and 
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TABLE I: Fourier coefficients from a fit of Eq. ((Tjl to tfie pion source radii squared. Tfie vaiues are in fm^. 



centraiity 


mean cent 


-Rout ,0 


-Rsidc.O 


o2 


-Rol.O 


^out,2 


^sido,2 


d2 

-'^long 2 


^os ,2 


0-2.5 % 


1.3 % 


29.34(10) 


24.33(08) 


33.10(11) 


-6.56(21) 


-0.25(25) 


-0.03(21) 


-1.43(26) 


0.19(17) 


Z.O-D /O 


O. ( /o 


zo.Ul(^Uo J 


Zo.Z\)[\JD ) 


ol.yi(^uy ) 


-0. 1 Z[io) 


-U.ODl^iy ) 


-U. 40(^10 ) 


-l.Zoyi)J ) 


A /I C /I Q\ 


5-7.5 % 


6.1 % 


26.62(08) 


22.04(06) 


30.61(10) 


-5.22(15) 


-0.96(19) 


-0.06(15) 


-0.59(19) 


0.46(12) 


7.5-10 % 


8.1 % 


25.39(17) 


21.00(13) 


29.22(21) 


-5.20(31) 


-1.15(32) 


-0.11(26) 


-1.25(32) 


0.63(21) 


10-15 % 


11.6 % 


23.83(19) 


19.87(14) 


27.31(23) 


-4.33(37) 


-1.14(35) 


0.42(27) 


-1.40(34) 


0.39(23) 


15-25 % 


17.5 % 


21.35(23) 


17.62(17) 


24.64(28) 


-3.67(48) 


-1.03(31) 


0.38(24) 


-1.27(30) 


0.35(21) 


25-70 % 


30.0 % 


14.28(49) 


12.29(35) 


17.42(61) 


-4.23(116) 


-0.98(30) 


0.36(23) 


-1.20(29) 


0.34(20) 



side directions indicate a pion source elongated out-of- 
plane. For comparison, the AGS [3 and RHIC [1, [11 
values, obtained by performing the fits using Eq. ([7]) on 
their published radii (the AGS results were subsequently 
corrected for their reaction plane resolution), are rep- 
resented in Fig. [1] by open symbols and stars. It ap- 
pears that the anisotropics in the out and side direc- 
tions at SPS energy are rather similar to those observed 
at RHIC. The geometrical pion source eccentricity can 
be quantified via e ~ 2i?gjjj, 2/-^side ^^'-^ ^ 
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FIG. 1: Azimuthal pion source eccentricity, represented by 
the second Fourier coefficient of squared radii 7?^($*), mea- 
sured in Pb-j-Au collisions as a function of centraiity. Positive 
and negative pion pairs have been combined. The mean pion 
transverse momentum is 0.23 GeV/c. For comparison, anal- 
ogous measurements at the AGS (open green symbols) and 
RHIC (blue diamonds and stars for 130 and 200 GeV, respec- 
tively) are shown. The last two panels show the first Fourier 
coefficients j and which should be zero by construc- 
tion. The estimated systematic error is 0.3 fm^. 



centrahty of 15-20%, a value of 0.043(27), significantly 
less than the initial fireball eccentricity Emitiai ~ 0.20 
The fact that the magnitude of the i?sidc anisotropy 
seems to be weaker than that of i?out indicates that not 
only the source geometry but also azimuthal dependence 
of the emission time might play a role. On the other 
hand, R'^^^ ^ ~ ^side 2+2-^os 2' averaged over ccntrahties, 
is 0.06(18) fm^, consistent with zero. This implies that 
the sum rule from [l^ . which is supposed to be valid for 
systems with emission time independent on the azimuthal 
angle, works rather well. 

The i?iong anisotropy is negative and roughly inde- 
pendent of centraiity. Series of checks were performed 
to make sure that this is not an artefact of the anal- 
ysis. This result might indicate that i?iong is sensitive 
to fluctuations of the azimuthal particle density. Hydro- 
dynamic calculation of central Fb-j-Au collisions at this 
energy yields a similar amount of anisotropy in -Rfo^g 2 
which, however, is centraiity dependent [l3|- This, and 
the fact that the hydrodynamic calculation ovcrpredicts 
the overall magnitude of i?iong, indicates that the knowl- 
edge about the mechanism leading to oscillations of i?iong 
is still incomplete. 

The source radius anisotropics for centraiity 15-20% 
are shown as a function of the collision energy in Fig. [2l 
The SPS result fits rather well into the beam energy sys- 
tematics. The apparent (albeit not statistically signifi- 
cant) fast change of R^^^^ 2 between AGS and SPS and, 
even more, the negative Rf^^^g ^ developing when going 

down in energy from RHIC to SPS, make the perspec- 
tive of a systematic study in the course of the low-energy 
scan at RHIC especially attractive. With the statistical 
errors of the present AGS data even significant structures 
in the energy dependence of the pion source anisotropy 
cannot be excluded. 



V. SUMMARY 

We have analyzed the azimuthal angle dependence of 
the pion HBT radii in Pb-j-Au collisions at the top SPS 
energy. The source anisotropy in the out and side di- 
rections has the same sign and similar magnitude as 
the one measured at the AGS and at RHIC, and in- 
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FIG. 2: Collision energy dependence of the pion source 
anisotropy in Au+Au and Pb+Au collisions at the 15-20% 
centrality. The meaning of the symbols is the same as in 
Fig. [T] The SPS results are rather similar to those obtained 
at RHIC, except for i?iongwhich is significantly off-zero. The 
estimated systematic error is 0.3 fm'^. 



dicates a pion source elongated out-of-plane. The side 
anisotropy is somewhat smaller than the other which sug- 
gests that finite emission times may play a role. The 
source anisotropy in the long-direction is negative for all 
centralities indicating that i?iong might be sensitive to 
particle density fluctuations. 
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